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VERTICAL STRUCTURE OF ADVECTION DOMINATED ACCRETION FLOWS 
Fateme Zahra ZeraatgarPand Shahram Abbassi 1,2 

ABSTRACT 

We solve the set of hydrodynamic (HD) equations for optically thin Advection Dominated Accretion 
Flows (ADAFs) by assuming radially self-similar in spherical coordinate system ( r,d,(j >). The disk is 
considered to be steady state and axi-symmetric. We define the boundary conditions at the pole and 
the equator of the disk and to avoid singularity at the rotation axis, the disk is taken to be symmetric 
with respect to this axis. Moreover, only the T rc f, component of viscous stress tensor is assumed and we 
have set vg = 0. The main purpose of this study is to investigate the variation of dynamical quantities 
of the flow in the vertical direction by finding an analytical solution. As a consequence, we found 
that the advection parameter, /“ , varies along the 9 direction and reaches to its maximum near the 

rotation axis. Our results also show that, in terms of no-outflow solution, thermal equilibrium still 
exists and consequently advection cooling can balance viscous heating. 
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1. INTRODUCTION 

Accretion onto a compact object such as black hole is a 
fundamental phenomenon in the universe and most likely 
is the primary power source in systems like X-ray binaries 
(XRBs), active galactic nuclei (AGNs) and gamma-ray 
bursts (GRBs). Observation of X-ray and gamma-ray 
emission lines from black hole accretion disks demon¬ 
strates the idea of forming a hot atmosphere above the 
accretion disk or perhaps radiatively inefficient flows (RI- 
AFs) (e.g. Blair et. al. 1984; Raymond 1993; Jimenez- 
Garate et. al. 2005.) 

In the standard accretion disk model (Shakura & Sun- 
yaev 1973) the energy released via viscosity is radiated lo¬ 
cally and the accreting flow becomes cool very efficiently. 
Therefore, this model cannot produce high energy spec¬ 
trum and the idea of existence of hot corona above disk is 
needed to predict such high-energy emissions. In terms 
of RIAFs, in fact, the generated heat through viscosity 
is stored as entropy and can be transported with flow 
inwardly rather than immediately being radiated away 
from the system. Consequently, the flow temperature 
becomes extremely high and nearly the virial tempera¬ 
ture. As a result, the disk can radiate high-energy emis¬ 
sion such as ganrma-ray (see Kato et al. 2008; Yuan & 
Narayan 2014 for more details). The disks and flows with 
this essential feature are named optically thin advection- 
dominated accretion flows (ADAFs). Historically, the 
significance of advection energy in hot accretion flows 
was firstly recognized by Ichimaru 1977 and more impor¬ 
tantly, a wide range of studies on optically thin ADAFs 
has been done by Narayan & Yi 1994, Narayan & Yi 
1995a,b and Abramowicz et al. 1995 where the disk was 
thermally stable (more details in Yuan & Narayan 2014). 

It should also be mentioned that in the field of hot 
accretion flows, a lot of numerical HD and MHD simula¬ 
tions have been carried out to investigate the dynamics 
of hot accretion flows and one of the most important 
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findings by those simulations is that the mass inflow rate 
decreases inwardly (e.g., Igumenshchev & Abramowicz 
1999, 2000; Stone, Pringle & Begelman 1999; Hawley, 
Balbus & Stone 2001; De Villiers, Hawley & Krolik 2003; 
Igumenshchev, Narayan & Abramowicz 2003; Yuan & 
Bu 2010; Pang et al. 2011; Yuan et al. 2012a,b; Bu 
et al. 2013) and it is not constant as considered before. 
Following those simulation results, one and two dimen¬ 
sional self-similar solutions of ADAfs in the presence of 
outflow and magnetic field have been done (e.g., Xu & 
Chen 1997; Blandford & Begelman 1999, 2004; Xue & 
Wang 2005; Akizuki & Fukue 2006; Abbassi et al. 2008; 
Zhang & Dai 2008; Bu et al. 2009; Jiao & Wu 2011, Ab¬ 
bassi & Mosallanezhad 2013; Mosallanezhad et al. 2013, 
2014, Samadi et al. 2015). Note that, the properties and 
dynamic of the hot accretion flows with magnetic field 
and outflow is beyond the scope of this work. 

In the optically thin ADAFs, the accretion rate is very 
low, M < 0.1 LEdd/c 1 2 where Lsdd is the Eddington lu¬ 
minosity and c is the speed of light. In addition, the 
optically thin ADAFs are geometrically thick disks, i.e., 
H/R < 1, where H is the disk’s height scale and R is the 
radius in cylindrical coordinates. It should be noted here 
that in the novel vertically average self-similar method¬ 
ology presented by Narayan & Yi 1994 (hereafter NY94), 
the energy equation expressed as, 

q adV =q + -q~= fq + (1) 

where q adv represents the rate of the entropy advection 
in the radial direction, i.e., pv r Tds/dR. Where, p is the 
density of the gas at the equatorial plane of the disk, 
v r is the radial velocity and also s and T are the spe¬ 
cific entropy and temperature of the gas respectively. 
In the equation ([1]), q + gives the total heat generated 
by viscosity per unit volume per unit time in the ra¬ 
dial direction. They defined the advection parameter as 
/ = q adv /q + which measures the fraction of the advec¬ 
tion energy stored as entropy. Consequently, (1 — /) will 
be radiated away from the system. They integrated the 
flow equation in the vertical direction. Making usual 
assumption such as steady state, axisymmety and a- 
viscosity they obtained a set of ordinary differential equa- 
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tions for the variables as a function of r. They have 
shown that the equations has a exact self-similar solution 
where all variables have power-law dependencies on r. 
Vertical integration is an standard approximation which 
has been used for thin disks where vertical thickness is 
usually much smaller than the local radius. While as we 
mentioned, optically thin ADAFs are geometrically thick 
disks, therefore, height-integrated approximation is not 
appropriate because the physical variables are not only 
function of r but also they should be function of the 
vertical direction, 9. Then, in the case of optically thin 
ADAFs, ID approach is not suitable. Later on, Narayan 
& Yi 1995a (hereafter NY95a), tried to solve the Hydro- 
dynamic equations (HD) in spherical polar coordinates 
(r, 9) with radially self-similar solutions but their results 
only corresponded to the simplest form of the advection 
parameter, viz., / = constant. 

To answer how the advection parameter varies along 
spherical polar angle, Gu et al. 2009 adopted a poly¬ 
tropic relation, p = Kp r , in the vertical direction which 
is normally used in vertically integrated geometrically 
thick disk model (e.g., Kato et al 2008). They defined 
the inclination 9 S near the polar axis regarding as the sur¬ 
face of the disk. Therefore, the main conclusion was that 
the optically thin ADAFs are geometrically thick since 
the free surface of the disk is very close to the polar axis. 
By taking into account the effect of a toroidal magnetic 
field and its corresponding heating, Samadi et al. 2013 
determined the thickness of advection-dominated accre¬ 
tion flows. Their results show that the vertical compo¬ 
nent of magnetic force acts in the opposite direction to 
gravity and compresses the disc; thus, compared with 
the non-magnetic case, in general the disc half-thickness, 
A 9, is significantly reduced. It should be emphasized 
that in both above mentioned works, the power index 
F considered with a typical value above unity, for ex¬ 
ample, r = 4/3. Also, the constant K is set equal to 
one, K = 1.0, to solve their one boundary differential 
equations starting from surface of the disk. Simulations 
carried out by De Villiers et al. 2005 revealed that the 
time average density drops faster than pressure in the 
vertical direction which means that the power index T 
should be less than one to satisfy the polytropic relation 
in the 9 direction. 

A few analytical solutions in the case of hot accretion 
have been presented and such solutions need to assume 
some simplifications. For instance, Shadmehri 2014 con¬ 
sider'd the energy equation of the gas with the inclusion of 
only T r( p component of viscosity in the viscous energy dis¬ 
sipation term and found a very creative analytical solu¬ 
tion. Later on, Based upon new simulation results men¬ 
tioned above, Gu 2015 (hereafter G15) repeated his pre¬ 
vious work, Gu 2009, with an original idea of finding an 
analytical solution. The behavior of physical variables 
are in satisfactory with those in NY95a except for the 
variation of isothermal sound speed and radial velocity 
profiles along the vertical direction. The main conclu¬ 
sion was, viscous heating and advection cooling cannot 
balance each other. Therefore, no thermal equilibrium 
exists under the purely inflow assumption. 

The analytical solution we present here is in the same 
methodology as described in G15, with three modifica¬ 
tions. Firstly, we define the first boundary conditions at 
the rotation axis 9 = 0° to increase the angular range of 


our calculation, we made this change because in terms 
of optically thin ADAFs, the disk is considered to be 
geometrically thick and there might exist low-dense with 
high temperature flow above the surface of the disk. Sec¬ 
ondly, following NY95a, to avoid singularity at the poles, 
the disk is taken to be symmetric with respect to this 
axis. The second change leads to find a relation between 
the value of constant I\ and the density at the rotation 
axis. Finally, we will adopt the modified ”a” descrip¬ 
tion of viscosity defined by Bisnovatyi-Kogan & Lovelace 
2007. In the next section we explain with more details 
why this form of viscosity is needed. With the aforemen¬ 
tioned modification, we can address whether there exist 
thermal equilibrium in the purely inflow case and check 
how the advection parameter change along the vertical 
direction. 

The outline of this paper is as follows. In section 2, 
the basic equations and boundary conditions are intro¬ 
duced. The numerical results are shown and discussed in 
more details in section 3. Finally, a brief summery and 
conclusions will be given in section 4. 

2. BASIC EQUATIONS AND BOUNDARY 
CONDITIONS 

2.1. Basic equations 

The standard hydrodynamic (HD) equations are em¬ 
ployed in spherical coordinate system (r, 9 , (p ) where the 
steady state accretion flow is taken to be axisymmet- 
ric (i.e., d/d(f> = 0). The gravitational potential of the 
central black hole is described in terms of Newtonian 
potential which is more convenient for the self-similar 
formalization, ip(r) = — ( GM)/r. In addition, the flow is 
in non self-gravitating regime and initially the relativis¬ 
tic effects are neglected. Following NY95a, we assume 
vg = 0 , which corresponds to a hydrostatic equilibrium 
in the vertical direction. Although, this assumption is 
not appropriate if you are investigating the effects of out¬ 
flow on the dynamics of accretion flow (see e.g., Jiao & 
wu 2011; Mosallanezhad et al. 2014 for more details). 
Therefore, the continuity equation and the three compo¬ 
nents of the equation of motion are as follows, 


^^^ 2pvr ) = 0 ’ 

( 2 ) 

- GM 1 dp 

r r 2 p dr' 

( 3 ) 

2 , 0 1 d P 

V * mte = - P w 

( 4 ) 

+ ^ = \ 1 (r 3 Tr+] ) ' 
r pr 6 or 

( 5 ) 


where v r and v^ are radial and azimuthal components of 
velocity, p is the mass density and p stands for the gas 
pressure. Besides, in equation <[5j), T r( j, represents the rep 
component of the anomalous stress tensor, respectively. 
It should be emphasized that in a real case, the mag¬ 
netic stress driven by the magneto-rotational instability 
(MRI) transfers the angular momentum outside the disk 
(Balbus & Hawley 1991, 1998). Since in our HD case 
we do not consider magnetic field, the anomalous shear 
stress tensor has been considered to mimic the magnetic 




Vertical structure of ADAFs 


3 


stress (see the HD simulations performed by Yuan et al. 
2012a for more details). This parameter can be written 
as, 

= (6) 

where p{= up) is the viscosity coefficient which deter¬ 
mines the magnitude of the stress and v is called the 
kinematic viscosity coefficient. There are a lot of un¬ 
certainties about how to prescribe such a viscosity pa¬ 
rameter. Most researchers adopt the “a” description for 
standard thin disks introduced by Shakura & Sunyaev 
1973 which is proportion to speed of sound as v = ahc s . 
Here, h is the disk’s height scale and a is a constant pa¬ 
rameter less than unity. We know that, if the viscosity 
coefficient scales with radius as v oc r 1 / 2 then the radial 
self-similarity will be possible. 

It should be point out that some simulations have been 
carried out for different forms of viscosity coefficients 
in accretion disks which found that the azimuthal com¬ 
ponents dominate other components (e.g. Stone et al. 
1996), therefore, in HD calculations, it would be more 
convenient to take into account the azimuthal compo¬ 
nents (see, Stone et al. 1999; Yuan et al. 2012a,b). 

NY95a obtained their solutions for ADAFs model cor¬ 
responding to the usual “a” description of viscosity. 
They also checked whether the results are sensitive to 
the viscosity by adopting different form as v = arc s and 
concluded that their solutions with new form are entirely 
similar to those with “a” description. This prescription 
may not be suitable in a real case, since there exists a 
low-dense corona above the disk with a nearly virial tem¬ 
perature. In addition, in the case of geometrically thick 
and hot disk (ADAFs model) the hottest temperature 
should be achieved at the rotation axis, 9 = 0,7r (see 
Figure 1 of NY95a for more details). On the other hand, 
the viscosity is due to MRI turbulence and this quan¬ 
tity should vanish at the surface of the disk (Bisnovatyi- 
Kogan & Lovelace 2007; Lovelace et al. 2009). Then, if 
the kinematic viscosity coefficient is proportional to the 
isothermal sound speed, v oc c s , then this quantity can¬ 
not vanish at the surface of the disk as MRI turbulence 
predicts. 

In order to avoid the disparateness in terms of the tur¬ 
bulent viscosity, following Lovelace et al. 2009, we adopt 
the modified ”a” description of viscosity as 

u = a^-g(9) (7) 

In the above equation, Qk (= y/GM/r 3 ) is the Keplerian 
angular velocity of the disk and g ( 9 ) is a dimensionless 
function which is equal to unity and zero in the body 
and surface of the disk, respectively (e.g., Lovelace et al. 
2009). For simplicity, we consider g (0) = sin# to satisfy 
the aforementioned conditions. 

We are interested in investigating whether the advec- 
tion parameter, /, which is normally considered to unity 
/ = 1, in the case of ADAF, remains constant along the 
polar angle or not (e.g., Narayan et al. 1995a; Xu & 
Chen 1997; Jiao & Wu 2011). Therefore, following Gu 
et al. 2009 and also Gul5, we apply polytropic relation, 
p = Kp T in the 9 direction as our last equation. Al¬ 
though, they obtained solutions by fixing K = 1 , we 


explain how this constant parameter will be determined 
in the next section. We also note that, the simulations 
of De Villiers et al. 2005 revealed that the power index 
r is less than unity. This is mainly important because 
their results show that the time averaged density drops 
faster than pressure from the equatorial plane to the po¬ 
lar axis. Based upon those results, T is set to be less 
than one throughout this paper. 

We adopt self-similar solutions in the radial direction 
to simplify the equations as 


P(r 1 9) = p(9)r 3/2 , (8) 


v r {r, 9) = \j^^-v r {9) = v K (r)v r (0), (9) 

v#(r, 6) = v K (r)v<p(9), (10) 

p(r, 9) = p(9)GMr~ 5 / 2 . (11) 

By substituting above self-similar solutions into Equa¬ 
tions 0-©, they will be reduced to 

- \vr{9? - V#) 2 = -1 + \Kp{ef-\ (12) 
u 0 (#) 2 cot# = KT P (9) r - 2 ^p-, (13) 

v r {9) = ~aKp(9f~ 1 sin#. (14) 


We will put Equations 03 and (fl4l) into Equation (fl3l) 
in order to obtain the differential equation for the den¬ 
sity. Before doing that, it should be note that, the first 
term on the left hand side of Equation m is very small 
compare to the other terms. This is because in the case 
of ADAF, radial velocity is very low and also the vis¬ 
cosity constant considered here is fixed as a = 0.1 (see 
equation (1141) 1. Hence, without any significant change in 
our results we can neglect this term and therefore, the 
differential equation will be written as (see G15 for more 
details): 


M°) 

d9 


cot# / p{9) 2 r 

~r~ ^ k 



(15) 


The above differential equation has an analytical solu¬ 
tion which will be obtained after introducing boundary 
conditions in the following section. 


2.2. BOUNDARY CONDITIONS 

We define the boundary conditions at the rotation axis 
# = 0, and equatorial plane, # = 7r/2 to occupy the 
angular range 0 < # < 7 t/ 2. Following NY95a, to avoid 
singularity at the polar axis, we assume the boundary 
condition as 

S = °: | = 0. (16) 


The above boundary condition leads to obtain the value 
of constant K. Therefore, in Equation (1151) . the term in¬ 
side the parentheses should become zero to satisfy above 
boundary condition. So, the K parameter can be derived 
as, 


F )Po r 


K = 


(17) 
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r 


per unit volume can be expressed as 

5 — 37 P V r 
2 (7 — 1 ) r ’ 


qddv _ 


+ __ 9a pv$ 


q = 


4 tvk 


sinO. 


(19) 

( 20 ) 


Therefore, by the vertical integration over q adv and q + , 
we can achieve Q adv and also Q + as 


Q adv = 2 f 2 q adv r sin 6 dd (21) 

Jo 

Q + = 2 I q + r sin 9 dO (22) 

Jo 

Then, the energy advection factor is given by f adv = 
Q adv /Q + . In the next section we will express the be¬ 
havior of all variables and also the variation of advection 
cooling and viscous heating by explanation and compar¬ 
ing our results to those in Gul5 and Gu 2009. 


Figure 1. Variation of K with T for different value of density at 
polar axis, p = 0.001, 0.01, 0.1. 


where po represents the value of density at the rotation 
axis. Fig. |T| shows variation of constant K versus T cor¬ 
responding to different values of po- Note that, since we 
will fix the value of the density at the equatorial plane to 
unity, p( 7 t/ 2 ) = 1.0, the magnitude of po in our Figures 
represents the ratio of the density at the polar axis to the 
mid-plane of the disk. It is clear from equation C3 and 
also Fig. |T| that K increases with increasing power index 
r for different fixed values of po- Another feature that 
can be seen in this figure is that the value of K cannot ex¬ 
ceed 2/5 which is the higher limit of this parameter with 
r index in a range 0 < T <1.0. As we explained before, 
we consider this range for T because according to the 
simulation results carried out by De Villiers et al. 2005 
the time average density drops faster than pressure in 
the 6 direction. This figure is compared with Figure (2) 
in G15 paper because at first, we decided to modify G15. 
Comparing Equation 03 with Equation (10) inside G15 
paper, obviously both will be equal if one considers the 
surface angle to be 0°. In order to avoid defining the pa¬ 
rameter A(= (po/poV'i'if) as the energy advection on the 
midplane of the disk (G15), we instead use the constant 
K which is the coefficient inside the polytropic equation. 

Now, we turn our attention to find analytical solution 
for the density in the vertical direction. By integrating 
Equation (fl5l) along the 9 direction, we can easily obtain 
the density profile as 


P(0) 


1 

5 K 


^ 5 /ip( 7 r/ 2 ) r 1 — 2 ^ sinfT 


■ 5(r~ij 



. ( 18 ) 

where the value of density at the equatorial plane is set 
to be p( 7 t/ 2 ) = 1.0 throughout this paper. To complete 
the specification of the results, we need to define the 
advection parameter, f adv . In the self-similar formalism, 
the advective cooling rate and the viscous heating rate 


3. NUMERICAL RESULTS AND DISCUSSION 

In this section, we first study the angular distribution 
of the dynamical quantities driven from Equations m- 
Q and Equation (fT51) as well. Fig. [2] shows the vari¬ 
ation of physical quantities with the polar angle, corre¬ 
sponding to T = 0.75, 0.85, 0.95 as three typical examples 
and the density at the polar axis as po = 0.01. The di¬ 
mensionless radial velocity, v t /vk , is plotted in the top 
left panel of Fig. [2] As it is seen, the radial velocity 
increases from the rotation axis to the equatorial plane. 
Actually, v r is zero at 9 = 0° and reaches to it’s maxi¬ 
mum at 9 = 7 t /2 for all T power index values. It should 
be note here that, in terms of F <1, G15 found that 
the radial velocity increases toward the rotation axis (see 
Fig. 1, panel (b) of G15). This contradiction is simply 
because, in this paper the modified form of a prescrip¬ 
tion for viscosity is adopted (see Equation (J7J and also 
Equation ([HI) for more details). In addition, our result 
shows that the larger T leads v r increases further. 

The top right panel displays the dimensionless az¬ 
imuthal velocity, i.e., v^/vk- It is clear that v$ increases 
from 9 = 0° to 9 = 7 r /2 and very close to the pole, v<f, 
vanishes and becomes zero. This behavior is clearly be¬ 
cause, as a second modification to G15, the disk is taken 
to be symmetric with respect to the polar axis. As a 
result, the centrifugal force will become zero at the rota¬ 
tion axis. Moreover, azimuthal component of velocity is 
larger for small values of T. In fact, the most variation of 
v<f, is belonged to the smaller value of T. As seen in this 
panel v$ is changed from v^ = 0 to just above v</, = 0.8 
for T = 0.75. 

The bottom left panel shows the vertical profile of the 
density p. It should be emphasized that the density pro¬ 
file is scale with the density value on the equatorial plane 
of the disk. In this figure, the minimum value of the den¬ 
sity at the rotation axis is considered as po = 0.01. As an 
overall trend, it is clear that the density increases from 
p ~ 0 at the rotation axis to p = 1 at the equatorial 
plane. What’s more, for the larger T as you can see, 
p increases sharper than the smaller one which means 
there exists extremely dense flow near the rotation axis 
and the disk is considered to be geometrically thick. 
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Figure 2. Variation of dimensionless physical quantities with polar angle 6. The dotted, dashed and solid lines correspond to T = 
0.75,0.85,0.95 respectively. Here a = 0.1 and po = 0.01. 


Finally, the bottom right panel of Fig. [2] shows the 
vertical variation of the isothermal sound speed, c 2 . As 
illustrated in this panel, c 2 has the decreasing trend from 
the rotation axis to the equatorial plane. Furthermore, 
from the figure it is clear that for the case with T = 
0.95, c 2 is almost independent of 8. Also the maximum 
variation of c 2 belongs to T = 0.75, from c 2 ~ 0.4 at 
8 = 0 (at nearly virial temperature)to about c 2 = 0.12 
at 8 = 7t/ 2. In this case, small value of T might be 
corresponding to the thin disk model with hot corona 
above the disk. Most importantly, our results are totally 
in satisfactory with those presented in NY95a. In fact, we 
should mention that by using an analytical solution for 
the case of no-wind self similar solutions, solving a system 
of complex ordinary deferential HD equations with two 
boundary condition it is not necessary (readers referred 
to NY95a for more details). 


As we explained before, the main purpose of this work 
is to check whether the advection parameter remains 
constant along the vertical direction or not. Fig. [3] 
represents the variation of the energy advection factor, 
f adv (= Q adv /Q +), with T on the range 0.5 < F < 1. In 
this figure the values of the density varies over the range , 
i.e., po = 0.001,0.01,0.1. It is shown that f adv increases 
with increasing T and reaches to the unity for high value 
of r corresponding to fully advection case. In contrary 
to G15 conclusions, this figure obviously demonstrates 
that in the case of no-outflow, i.e., ve = 0, there exists 
thermal equilibrium and therefore, advection cooling can 
balance viscous heating. Also as it is seen, for three typ¬ 
ical values of the density on the rotation axis, the fully 
advection take place when F ~ 0.95. Therefore Fig. 1 
together with this figure show when the energy equa¬ 
tion is replaced with the polytropic equation of state, 
p = Kp T , the fully advection case will be possible if con- 
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Figure 3. Variation of the energy advection factor, f adv (= 
Q adv /Q+), with r. The dotted, dashed and solid lines correspond 
to po = 0.001, 0.01, 0.1 respectively with a = 0.1. 



Figure 4. Vertical profile of the energy advection factor, /(0)(= 
qadv /q+y The dotted, dashed and solid lines correspond to po = 
0.001,0.01,0.1 respectively. Here, T = 0.95 and a = 0.1. 


stants K and F vary only in the range 0 < K < 0.4 and 
0.5 < T < 1 respectively. 

Finally, the vertical variation of advection factor, 
f(9) = q adv /q ms , (the ratio of advection cooling to the 
viscose heating per unit volume) is plotted in Fig. [I] 
for the same three values of the density in Fig. [3] with 
F = 0.95. It can be seen that the value of the energy 
advection factor is not constant in the 9 direction and 
increase from equatorial plane towards the rotation axis. 
In addition, the lower value of the polar axis density 
causes to the smaller value of f(9). It is also seen that 


the energy advection factor is below unity, i.e., f{9) < 1 
for p = 0.001. What’s more, for the case with p = 0.01 
the value of the energy advection factor is first below 
unity an then become greater than one near the rotation 
axis. Furthermore, the larger value of the density close 
to the rotation axis, i.e., p = 0.1, causes to the upper 
limit of f(9), i.e., f(9) ~ 2.8 at 9 = 0 and more impor¬ 
tantly this factor is always above unity for all vertical 
range. Therefore unlike NY95a assumption which has 
been considered / = 1.0 throughout the angular direc¬ 
tion, our result clearly indicate that the advection pa¬ 
rameter is a function of 9 , i.e., f(9), and may exceed 
unity in some cases. 

4. SUMMERY AND CONCLUSIONS 

In this paper, we tried to solve the hydrodynamic 
equations of optically thin ADAFs in spherical coordi¬ 
nate system (r, 9 , (f >) where the steady accretion disk is 
considered to be symmetric with respect to the rota¬ 
tion axis as well as the equatorial plane. The central 
black hole gravity is described as the Newtonian poten¬ 
tial, since this form is more convenient in the self-similar 
solutions, tl>(r) = — ( GM)/r. In addition, instead of us¬ 
ing energy equation with constant value for the advec¬ 
tion parameter, /, following Gu et al. 2009 and G15, we 
adopted the polytropic relation in the vertical direction 
as p = Kp T . Compared with G15, we made three mod¬ 
ifications. Firstly, the vertical range of the calculation 
is enhanced from rotation axis 9 = 0° to the equatorial 
plane of the disk, 9 = n /2. This change have been made 
since in optically thin ADAFs, the disk is geometrically 
thick, i.e., H/R < 1. Secondly, following Narayan & Yi 
1995a to avoid singularity at the pole, the disk is taken 
to be symmetric with respect to this rotation axis. This 
change causes to find a relation between the value of 
constant I\ and the density at the rotation pole. Finally, 
the modified ”a” description of viscosity is adopted, (see, 
e.g., Lovelace et al. 2009 for more details). 

By the above mentioned modifications and following 
G15 methodology, we could find an analytical solution 
for optically thin ADAfs. The presented results showed 
that unlike G15, the radial velocity decreases towards 
the rotation axis for all T values. In addition, be¬ 
comes zero at rotation pole and this is due to the second 
modification. Furthermore, c^, the same for NY95a, has 
a decreasing trend from 0 = 0 (virial temperature) to 
9 = tt/2. 

Besides in contrast to G15, our solution represents 
the existence of thermal equilibrium in vertical direc¬ 
tion without outflow emanating. So advecting cooling 
can balance viscous heating effectively. Moreover, If a 
polytropic relation is used rather than the energy equa¬ 
tion in the vertical direction for the fully advection, K 
and F should be only on the range 0 < K < 0.4 and 
0.5 < r < 1, respectively. At last, the value of the en¬ 
ergy advection factor is not constant in 9 direction and 
increase from equatorial plane towards the rotation axis. 

In spite of the simplicity of our model in viscosity 
and the disk itself, we think that the presented semi- 
analytical results give us a better understanding of such 
a complicated system. It is good to note here that some 
modifications can be applied to ameliorate this study. 
Regarding to the radial self-similar approximation, it can 
not ensure us that this solution is definitely relevant to 
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the real accretion flows. In addition, the Newtonian po¬ 
tential was taking into account rather than Paczynsky 
& Wiita potential to avoid the general relativity effects 
in the innermost region of accretion disk. Not only the 
r<p component of the viscous stress tensor, T r< j,, should 
be employed, but also the other components such as Tg^ 
ineluctably should be taken into consideration. How¬ 
ever, the advection factor f adv in the energy equation 
was found to be function of the vertical direction, but 
this parameter must be a function of radial and more sig¬ 
nificantly the mass accretion rate of the disk, i.e., m. To 
resolve the aformentioned remarks, considering vg and 
also component of stress tensor may cause to have 
more promising results in the future works. 

The authors would like to thank Amin Mosallanezhad 
for his useful suggestions and discussions. We also would 
like to appreciate the referee for his/her thoughtful and 
constructive comments in the early version of the pa¬ 
per. This work was supported by Ferdowsi University of 
Mashhad under the grant 3/37875 (1394/04/03). 
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